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ABSTRACT 



This study describes the continuing development of a computer model for the 
simulation of burst kill probabilities for air defense gun systems firing projectiles at 
maneuvering aircraft. The computer simulation developed by Keeling is modified by 
adding a high explosive proximity and contact fuzed round to the simulation. The 
objective of this thesis is to develop and analyze the shipboard anti-air defense problem 
in order to choose the best air defense gun system. The air defense gun system studied 
consists of a 40 mm gun, firing fragmenting proximity and contact fuzed projectiles, an 
early warning radar system, and a fire control system. The aircraft vulnerability and 
the gun data used in this thesis are entirely synthetic to avoid security classification 
problems. 
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I. INTRODUCTION 



A. BACKGROUND 

In modern times it often occurs that during battle, friendly ships are threatened 
by enemy aircraft. This thesis focuses on a ship's ability to defend itself against aircraft 
with a 40 mm gun system. The computer program presented in this paper is a 
modification of a simulation written in the A Programming Language(APL) by Capt. 
C. Keeling. [Ref. 1] Keeling’s program simulates engagements between air defense gun 
systems firing non-explosive projectiles at maneuvering aircraft. In his thesis, Keeling 
suggested further research on proximity fuzed and contact fuzed high explosive(HE) 
rounds. This paper expands Keeling's simulation to include the HE round with both a 
contact fuze and a proximity fuze 

The simulation shows how changing gun parameters can affect a gun system's 
probability of kill against an aircraft. These parameters include: target range, projectile 
size, ballistic dispersion, muzzle velocity, firing rate, radar error, linear fire control, type 
of fuze (a new feature) and aircraft probability of kill given a hit or warhead lethal 
radius (a new feature). The use of this program increases the information available on 
the effectiveness of gun systems with and without HE warheads. The ultimate 
measure of success in the design and performance of an air defense gun system is the 
ability of the system to shoot down low-flying hostile aircraft within the combat 
environment. 

B. OBJECTIVE 

The objective of this thesis is to continue the development of a computer model 
that will be helpful in making decisions regarding the design and operation of an air 
defense gun system. In this thesis, the program is used to gain insight into how each 
of the design parameters of the gun system affect the probability of kill and which of 
the parameters are most important in the gun system design for air defense. In 
particular, the influence of the HE round on gun effectiveness will be studied. 
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C. APPROACHING THE PROBLEM 

The gun system is located on a ship at the origin of an X, Y, and Z coordinate 
system. The playing area for the model is the quarter of a hemi-sphere entirely within 
the first quadrant of the X, Y and Z plane as shown in Figure 1.1. The size of the 
playing area may be adjusted by modifying the program's source code. 




Figure 1.1 Playing Area of The X, Y & Z Plane. 
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The model assumes one-on-one engagements, no suppressive fire from the aircraft, 
and no terrain effects. The target is not necessarily attacking the gun system and 
may or may not maneuver while within the gun's effective range. 



D. ORGANIZATION OF THE MODEL 

The air defense system consists of an early warning system, a tracking radar, a 
weapon control system and a gun with one or more barrels. The model studied here is 
a twin 40 mm gun using high explosive rounds located onboard the ship at the origin 
of the playing area. As an aircraft enters into the radar detection envelope of the early 
warning radar, the aircraft is detected and identified as hostile by the ship. The 
tracking radar is given the target's position. This radar tracks the aircraft and 
continuously computes the target range, elevation, azimuth, and velocity. The gun fire 
control computer solves for the lead angle in azimuth and elevation based on the 
calculated target positional data assuming no target maneuvering. The gun is 
positioned at the lead angle and fires a burst of proximity or contact fuzed high 
explosive rounds. The projectiles fly out toward the moving target. The target 
continues on its randomly selected flight path. The miss distance and assumed ballistic 
dispersion are used in the model to determine the probability of aircraft kill for each 
round Pj^sS" ^ or P rox * m i t y fuzed warheads, the projectiles do not have to hit the 
target to cause damage. The probability of aircraft kill given a detonation, P^ ^,is 
determined using the Carlton diffused Gaussian kill function . For the contact fuzed 
warheads, the projectiles must hit the target in order to cause damage. The probability 
the aircraft is killed given the hit on the target, Pj^ is represented by the two- 
dimensional Carlton hit function. 

1. Scenario 

The general approach in this model is the development of the kill probability 
for one gun system firing one burst at one randomly maneuvering target. The 
probability of kill for each round is used to compute a burst kill probability. The gun 
system computes target flight path characteristics, predicts a time-dependent intercept 
point in space, positions the gun and fires the projectiles. For contact fuzed HE 
rounds, the target vulnerability is represented by a vulnerable area (A v ) in a plane in 
space that is perpendicular to the slant range (R) between the gun and the center of 
the target. For the proximity fuzed warhead, the target vulnerability is represented by 
a lethal radius. 
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2. Model Input 

Inputs are required for the gun, the fuze, the target, and ranges of 
engagement. Gun inputs are muzzle velocity, projectile coefficient of drag, burst size, 
rate of fire, and angular ballistic dispersion. The fuze input is the probability of fuzing, 
which is independent of range. Aircraft inputs are velocity, altitude, acceleration 
performance capabilities, and vulnerability data based on the type and size of 
projectile and fuze being used. A flight path generates, for several different types of 
aircraft, maneuvers using a Monte Carlo simulation. 
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II. PROXIMITY FUZED WARHEAD 



A. BASIC THEORY 

A proximity fuze detects the target and causes detonation of the high explosive 
warhead in the vicinity of the target. A warhead usually contains a powerful but 
relatively insensitive high-explosive that can only be initiated by the heat and energy 
from the primary explosive in the fuze [Ref. 3]. Proximity fuzes accomplish their 
purpose through "influence sensing" with no contact between the projectile and the 
target . These fuzes are actuated by some characteristic feature of the target, such as a 
reflected radio signal, an induced magnetic field, an interrupted light wave, a pressure 
measurement, or an acoustical impulse. "Proximity sensing" results in detonation of 
the bursting charge in the vicinity of the target or targets. A direct hit is not necessary 
to disable an aircraft and achieve the desired effect. An aircraft may be damaged in 
any part of the large volume occupied by fragments from the exploding round. Since 
several air targets may exist simultaneously, the sensing equipment must be able to 
isolate a selected target. Proximity fuzing had its origins in England early in World 
War 11. Proximity fuzing is represented as VT fuzing (a code name used during Word 
War 11 to imply variable time fuzing) [Ref. 2 p. 76]. 

B. TARGET SUSCEPTIBILITY 

The modeling and quantification of the individual events and elements in an air 
defense encounter is referred to as a susceptibility assessment. Susceptibility is the 
inability of an aircraft to avoid being damaged in the pursuit of its mission, to its 
probability of being hit [Ref. 2] The level of susceptibility of an aircraft in encounter 
with a threat is dependent upon three major factors: 

• The scenario includes the physical environment in which the encounter takes 
place, the air defense deployment and activity, and the aircraft flight path and 
tactics, including any supporting forces. 

• The air defense characteristics, operations, and lethality. 

• The aircraft observables or detectable signatures, any countermeasures used, the 
aircraft performance capabilities, and any self-protection armament are the 
important factors associated with the aircraft itself. 
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One of the most important susceptibility measures is the closest point of 
approach or miss distance of a projectile to a maneuvering target. 

In general, the smaller the miss distance, the more likely the aircraft will be hit. 
Whether or not the aircraft is hit depends upon the presented area of the aircraft. This 
study assumed the presented area Ap of the aircraft is based on the six cardinal sides. 
The most the projectile can see in the intercept plane is part of three sides. The threat 
aircraft used in the study has a presented area from the area of the top A t0 p (50 m , 
equal to the bottom), the area of the side A^g (30 m the same on each side), and 
the area of the front Af ront (5 m^, equal to the rear). 

C. MISS DISTANCE 

As shown in Figure 2.1, the CPA is the mean value for an assumed Gaussian 
distribution of the projectile miss distance and the ballistic dispersion is the round to 
round variance about that point. 




Figure 2.1 Ballistic Dispersion, Bias, and CPA. 



The miss distance form for the round is given by 
M bx 2 + % 2 = CPa2 (eqn 2.1) 
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where 

• ^b* 'l*b v ) is the mean or bias in the x and v direction. 

D. THE FUZE 

The function of the fuze is to initiate the detonation of a warhead at a time and 
place such that maximum damage will be inflicted upon the target. Since the warhead 
damage volume is preset, the achievement of maximum damage depends on the time of 
fuze initiation (Ref. 3]. The fuze package consists of a safety-and-arming device to 
keep the weapon safe until it is deployed and clear of friendly forces, a detonator to 
initiate the HE charge detonation, a device that senses the presence of a target (known 
as the target detection device (TDD)), and a logic circuit that initiates detonation at 
the proper time. 

The functioning of a fuze can be represented by the probability P f which may or 
may not be dependent upon the miss distance. For example, a proximity fuze with an 
S0% probability of detonating the warhead within a miss distance range r of 100 ft 
and a 0% probability beyond r, can be modeled in the form 

Pjfx.y) = 0.S when(x“+ y“)^’“ < 100ft 
Pjfx.y) = 0 when (x~ + v^)^“ > 100ft 

E. TARGET VULNERABILITY 

Vulnerability refers to the inability of the aircraft to withstand one or more hits 
by the damage mechanisms, to its vincibility, to its liability to serious damage or 
destruction when hit by gun fire [Ref. 2]. The critical components on an aircraft are 
those components which, if either damaged or destroyed, would lead to an aircraft kill. 
A general procedure has been developed for determining the critical components, their 
possible damage or failure modes, and the effects of the component damage or failure 
upon the continued operation of the aircraft. The procedure consists of 

• A selection of the aircraft kill levels or categories to be considered. 

• An assembly of the technical and functional description of the aircraft. 

• The determination of the aircraft components of the aircraft and their damage 
caused failure modes for the selected kill levels. 
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1. State Transition Matrix Method 

The study by Keeling used the state transition matrix method to describe the 
damage state of the aircraft after one or more hits. This method allows the 
consideration of redundancy in the critical components and accounts for the fact that a 
single shot cannot kill more than one redundant critical component. This method 
states that a sequence of independent events can be modeled as a .Markov process. 
The state space of the Markov process includes kill of the aircraft, no kill of the 
aircraft, a state space for a kill of each of the critical redundant components, and 
combinations of kills of redundant critical components that eventually lead to a kill of 
the aircraft., The state of the target at the begining of the engagement is represented 
as the initial state vector S f) . One transition represents the state of the target after one 
random hit on the aircraft, and + * represents the state of the target after the hit. 
The matrix that accounts for the transition is |TM|. Thus 

S i +I = !TM|si (eqn 2.2) 

The state transition matrix method is not appropriate for the HE round. 

2. Vulnerability Model for Externally Detonating Warheads 

For the proximity fuzed HE rounds, the aircraft's vulnerability can be 
represented by a q function of x and v proposed by Carlton 

P K D = ex P(-* 2 r 0 2 )*exp(-y 2 /r 0 2 ) (eqn 2.3) 

The scaling parameter Tq can be related to the warhead lethal radius q, such that when 
r=r h P K D = ^' Collocating the one-dimensional Carlton kill function through 0.5 
when r= ^ gives Tq= 1.20 r| as shown in Figure 2.2. 

F * P KSS 

This simulation uses the circular Carlton kill function method as presented by 
Ball (Ref. 2 p. 320]. The equation for P^g is 

P KSS = V P f( al ^ /2 )*exp(-0 l bx 2 “H^by 2 'P» ( ec l n 2 - 4 ) 

where 
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Figure 2.2 Function of X and Y. 

• a is ( 2<r dx 2 + r 0 2 ) 

• P is (2cr dy 2 + r Q 2 ) 

• (cr dv <T dv ) is the ballistic dispersion of the round in the x and y direction. 

• r 0 2 is a scaling parameter of the warhead lethal radius r t such that r 0 = 1.2r.. 

• P f is the probability of fuzing (assumed constant). 

The calculation of the probability of kill for a variety of ranges is made simpler 
by using angular errors for the dispersion errors. From the Figure 2.3 it can be seen 
that a x j and a x 2 are the angles representing the angular standard deviation of the 
dispersion errors about the closest point of approach. Figure 2.3 shows that 

^dx = a dx R ^ n 2>5 > 

where 

• a dxl = a dx2 = a dx 

• otj x is the angular ballistic dispersion in the x direction. 

By the same theory in the y direction can be shown to be 



IS 



^dy a dy^ 

where a dy is the angular ballistic dispersion in the y direction. 




Figure 2.3 Effect of Angular Dispersion with Range. 



Furthermore assumption is made that the angular ballistic dispersion as shown in 
Figure 2.2 is the same in both direction, hence 

«dx = a dy = « 

The equation for the Carlton kill function model. Equation 2.3, simplifies to 

P KSS = < r o 2f Y V) * exp(-(H bx 2 + M by 2 )/(Y)) (eqn 2.6) 

where 

• y is (2 a 2 R 2 + r Q 2 ) 

The final form for the probability of kill is obtained by using Equation 2.1, giving 
P KSS = ( r o 2l V Y) ,1 ‘ e xp(*CPA 2 / Y) (eqn 2.7) 



For this study, rj is taken as five meters. Thus, Tq is six meters. 
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G. PROBABILITY OF KILL GIVEN A BURST 

The probability of an aircraft kill for a burst of rounds is directly proportional to 
the probability there is a clear line of sight from the detecting element to the aircraft. 
Given that the aircraft has been detected, the probability that one or more propagators 
will be fired at the aircraft is detonated by Pp. Thus, the probability that the aircraft is 
killed by one shot in an encounter, Pj^/p is given [Ref. 2. p. 320] by 



• Pq is the probability that the aircraft has been detected from the start of a 
search up to the present time t. 

• Pp is the probability that a propagator will be fired at the aircraft. 

• This study assumes that Ppj* Pp = 1. 

The probability that the aircraft is killed by the sequence of N gun shots is unity 
minus the product of the individual probabilities of survival for each shot. Thus, for 
the burst of N rounds P^ss in Equation 2.8 is replaced with 




(eqn 2.8) 



N 




(eqn 2.9) 



»=1 

where the subscript i denotes the ith shot. 
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III. CONTACT FUZED WARHEAD 



Contact fuzes detonate the HE warhead upon contact with the target. 
Detonation can be made to occur on the surface or inside the target, depending upon 
the design of the fuze. 

A. TARGET SUSCEPTIBILITY 

The presented area A p of the aircraft is defined to be the area of the aircraft seen 
by the projectile in a plane normal to the path of the projectile at intercept. The 
procedure for determining the presented area seen by the round is determined based 
upon the difference between the projectile's velocity and the relative velocity of the 
projectile with respect to the aircraft. This procedure is described in detail in [Ref. 1]. 

B. TARGET VULNERABILITY 

For warheads that do not use a proximity fuze, such as small arms and HE 
warheads with contact fuzes, a hit on the aircraft must occur to cause damage. Thus, 
the kill function that defines the probability the target is killed due to a propagator 
becomes ^ for random hits, and integration is carried out over the extent of the 
aircraft. In this approach, the probability the aircraft is hit by the propagator, Ppj, is 
computed and multiplied by the probability the aircraft is killed given the random hit 
on the aircraft, Pjt pj. The computation of Pj_j is developed by Keeling and 
accomplished using the Carlton hit function. The aircraft vulnerability is represented 
by its vulnerable area Ay, centered at the aim point, and any hit on the vulnerable 
area causes a kill. The computation for the vulnerable area is the product of presented 
area Ap and the probability the aircraft is killed given the random hit on the aircraft 

Pk.H l Ref * 21 

Ay= H*Ap ( ec l n 3-1) 

This model assumes the Pj^ pj value for the round is 0.25, and Ay is given as 
Ay= 0.25* Ap (eqn 3.2) 
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C - P KSS 

This simulation uses the Canon kill function for the P^^ as presented in [Ref. 

2 ]- 



p - 
r KSS 



2 2 1/2 2 2 1/2 
(2no x +x 0 ) (2 kg +y Q ) 



exp 






(2kg r 



■*o> 






(2KG* + y 2 0 )) 



(eqn 3.3) 



The product of x Q and y Q can be taken as (l)equal to the presented area, Ap, or (2) 
equal to the vulnerable area, Ay (see Ref.2, PP 316-315). This study assumes Xq = \'q 
and uses Ap for the product. The means and deviations in Equation 3.3 are 
determined in the same manners as was done for the proximity fuzed round in Chapter 
II. Thus, the final form for the probability of kill given a single shot becomes 

P KSS = (.25)Ap,(27ta 2 R 2 + Ap 2 )*exp{-7rCPA 2 /(2na 2 R 2 + Ap 2 )} (eqn 3.4) 



D. PROBABILITY' OF KILL GIVEN A BURST 

When multiple shots are considered, the single hit Pj^/jq is used for each shot. 
To determine the probability of an aircraft kill for a burst of rounds, the same 
procedure uses for the proximity fuzed warhead is used. Thus, 

N 

/Burst = * ~ I! ( * - ^ eC l n 

1=1 
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IV. COMPUTER PROGRAM 



Keeling's program is modified by adding subroutines PROXFUZE, CONTFUZE, 
BIRSTKILL, and by changing the appropriate data to a 40 mm gun. This computer 
program consists of a main program and utilizes fifteen subroutines. The program was 
written in the A Programing Language (APL). The computer program listing is shown 
in APPENDIX A, and a flowchart of the program is depicted in Figure 4.1. 



A. MAIN PROGRAM 

The main program uses a DRIVER program for initiating simulation and probes 
the user for inputs. The program DRIVER calls the GUN subroutine for simulating 
gun fire solution. After the simulation has been initiated, subroutine GUN simulates 
the air defense engagement by utilizing the appropriate subroutines necessary to 
simulate the engagement of an aircraft by the on-board ship air defense gun. 
Subroutine GUN collects statistics on the probability of kill and where the actual slant 
range of intercept occurred. These data are then sent back to the main program 
DRIVER. 

B. SUBROUTINES 

1. Initiation 

Subroutine GUN simulates the air defense engagement by calling the 
appropriate subroutines necessary to simulate the engagement of the aircraft by the on- 
board ship air defense gun. 

2. Generation 

The generation of a situation for air defense, aircraft's position, the radar 
located position with radar errors, and flight path maneuver of aircraft is subroutines 
RANGECHEK, RADAR, and AIRCRAFT. 

3. Calculation 

Subroutines INTERCEPT, MISSDIST, and MISSVECTOR are used to 
determine the predicted intercept point for the median round in the burst and the 
closest point of approach for each round in the burst. 
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Figure 4.1 Flowchart of Computer Program.. 



4. Susceptibility 

Subroutines APDATA, PKVECTOR, CONTFUZE, and PROXFL’ZE are 
used in determining the susceptibility of the aircraft. APDATA determines the 
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presented area of the aircraft that is used in the Carlton kill function model in 
subroutine CONTFUZE. PKVECTOR stores the data for each round in the burst. 

5. Vulnerability 

Subroutines PKVECTOR. CONTFUZE. PROXFUZE. and BURSTKILL are 
used in determining the vulnerability of the aircraft. Subroutine PKVECTOR stores the 
data for each round in the probability of kill. Subroutine CONTFUZE and 
PROXFUZE are used in determining the probability of kill for a single shot whether or 
not the projectile hit the target. 

6. Random Number 

Subroutines U'NTF, UNTRAN'D, and NORRAND are used to generate the 
random number required in the simulation. Subroutine UNTF generates uniform 
random numbers in the interval(a, b). UNTRAN'D generates pseudo-uniform random 
numbers. NORRAND is the generation of normal random numbers. 

C. SUBROUTINE PROXFUZE 

Subroutine PROXFUZE is used in determining the probability of kill for a single 
shot of the HE proximity fuzed round. This vulnerability data is used by subroutine 
BURSTKILL to determine the probability of kill given burst. 

D. SUBROUTINE CONTFUZE 

Subroutine CONTFUZE is used in determining the probability of kill for a single shot 
of the HE contact fuzed round. This vulnerability data is used in subroutine 
BL'RSTKILL to determine the probability of kill given burst. 

E. SUBROUTINE BURSTKILL 

Subroutine BURSTKILL determined the probability of kill for a given burst. 
This probability of kill data uses the vulnerability of the aircraft. 



25 



V. ANALYSIS AND CONCLUSIONS 



This program model was initiated with the realization that the probability of kill 
is one measure of effectiveness that has a broad base application in the field of air 
defense gun system evaluation. This program for determining gun system probability 
of kill has simulated several parameters. These parameters describe the gun, target 
environments, and engagement procedures whereby one gun system encounters one 
maneuvering target. The data used is given in APPENDIX D. This study uses a box 
plot to illustrate the effects of several important parameters. This graphical 
imformation gives a quick sense for the distribution of Pj^ for the individual test 
values of each of the parameters. The vertical box contains the middle fifty percent of 
the distribution, and the individual dots show all the outliers. The mean value of the 
probability of kill at each of the test value settings are connected with a line. The result 
of each simulation is listed with the Empirical Statistics table in APPENDIX E. 

A. TYPE OF FUZE 

The probability of kill of the aircraft in a one-on-one encounter with both the 
proximity fuzed warhead and the contact fuzed warhead, P^ £. is shown in Figure 
5.1. The effects that a proximity and a contact fuzed warhead have on the lethality of 
the 40 mm gun using a HE rounds can be seen in this Figure. The proximity fuzed 
warhead shows a higher probability of kill given a burst than the contact fuzed 
warhead against a maneuvering target. Consequently, the remaining studies 
concentrate on the proximity fuzed warhead. 

B. BALLISTIC DISPERSION 

Figure 5.2 shows an example of how the 40 mm gun probability of kill changes 
with angular dispersion. The degree of dispersion can immediately be seen, as well as 
how many of the data points were considered as outliers. Since the dispersion settings 
for optimality is arbitrary in this range, the assumption is made that the dispersion 
angle for the gun is 1.0 milliradians. 

C. RANGE EVALUATION 

The evaluation of the Pj^ for different range values was conducted for the 
maneuvering target. The maneuvering target was engaged by one burst for each range 
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lethality or proximity fuze 




LETHALITY Of contact fuze 




Figure 5.1 Graphical Comparison of Fuze on 40 MM gun Lethality. 

value. The gun produced the highest probability of kill at the 700 meters, and the 
effective range was 1300meters against the maneuvering target. The air defense for a 
real situation requires about 2000meters for successful engagement. 
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EFFECT OF AWCULAK WSPCRSION 




Figure 5.2 Effects of Ballistic Dispersion on 40 MM gun lethality. 

D. EFFECTS OF TARGET MANEUVERS 

Figure 5.4 shows the effects that a maneuvering and non -maneuvering target 
have on the lethality of the 40 mm gun. Against the non-maneuvering aircraft, the 
probability of kill given a burst is about twice that for the maneuvering target. That 
means that the lethality of the air defense gun requires a non linear fire control 
predictor based upon the tactical air combat maneuver(TAC.M). 

E. EFFECTS OF RADAR ERROR 

Figure 5.5 shows the effect of radar error on the 40 mm gun. In today's combat 
environment, the radar accuracy is an important factor. That means a large 
improvement in lethality can be obtained by reducing the radar error. The 100 percent 
error represents the current state of the radar system, so greater or less than 100 
percent reflects great differences in lethality due to technological improvements to 
today's radar system. 
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EFFECTS OF RAMCES ON 40 MM GUN 




Figure 5.3 Effects of Range on 40 \1M Gun Lethality. 

F. EFFECTS OF MUZZLE VELOCITY AND RATE OF FIRE 

Increased rate of fire and muzzle velocity can significantly improve the 
probability of kill for the rounds in the burst, as shown in Figures 5.6 and 5.7. 
Increasing the rate of fire decreases the time between projectiles, which means there are 
more projectiles in the area of where the aircraft is at intercept. Thus, the probability 
of kill is increased for each round in the burst. Increasing muzzle velocity decreases 
the time required to get to the target, and thereby reduces the effects of target 
maneuver. 

G. CONCLUSIONS 

This study has been conducted from a theoretical and intuitive point of view with 
conclusions withheld in many cases pending validation. The execution of this program 
demonstrated what appear to be some clear improvements is gun lethality. The sample 
results in each of the graphs show the difference in the probability of kill values. They 
represent a convenient and practical base from which to study or perform parametric 
analyses on existing systems or on systems still on the drawing board. The 40 mm gun 
with the proximity fuzed warhead considered here for the given scenario and target 
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effect of maneuvering target 




EFFECT OF NON-MANEUVERING TARGET 




Figure 5.4 EfTects of Maneuver and Non-maneuver on 40 MM Gun lethality. 

vulnerability data has an effective range of about 1300 meters against a maneuvering 
target. This gun can be improved by improving the radar, the fire control computer, 
projectile muzzle velocity, gun system rate of fire. These improvements should be 
attempted in the order given. 
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EFFECT OF RADAR ERROR 




Figure 5.5 Effects of Radar Error on 40 M\1 Gun Lethality. 

Another air defense gun system study should be conducted to determine the 
synergistic effects of changing more than one parameter simultaneously. The model 
should be further developed to account for non linear fire control. This can be easily 
done by changing the radar and intercept section of this program. 
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EFFECT Of MUZZLE VELOCITY 




Figure 5.6 Effects of Muzzle Velocity on 40 MM Gun Lethality. 
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EFFECT OF FIRING RATE 




Figure 5.7 Effects of Firing Rate on 40 MM Gun Lethality. 



APPENDIX A 

APL COMPUTER PROGRAM 
1. KEELING'S PROGRAM 
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[17] 
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[36] 
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[45: 

[46] 



V DRIVER-, COUNT', I 

q************************************************** 
a******* DATA INPUT ******************** 

q ************************************************** 

A ' HOW MANY DATA POINTS DO YOU WANT PER RANGE? ' 

DATA* 50 
DATA+DATA-1 

a ' WHAT IS THE CUN FIRING RATE IN ROUNDS PER MINUTE? ' 

FI RE RATE*- 1000 

a 'WHAT IS THE MUZZLE VELOCITY? ' 

70^1150 

a ' WHAT IS THE BURST SIZE? ' 

BURST*- 5 

a ' WHAT IS THE ANGULAR DISPERSION IN RADIANS ? 1 
DISP* 0.0015 

a 'WHAT IS THE DRAG FACTOR B = CDxV0*0 . 5 ? ' 

B*0 . 15 

a ' WHAT ARE THE RANGES YOU WANT TO LOOK AT? ' 

RANGES * 1900 1300 700 100 
a 1 WHAT ARE THE AREAS OF THE AIRCRAFT? ' 

AFRT* 5 
ASIDE* 30 
ATOP* 50 

a * * * PROGRAM EXECUTION AND DATA COLLECTION ******* 
q*** P IS PROBABILITY OF KILL ****** 

A * * * R SLANT IS RANGE OF TARGET AT INTERCEPT * * * 

q*************************************************** 

R*l 0 
P«-l0 
1*1 

CO UNT*p RANGES 
LOOP : RANGE*RANCES [I] 

GUN 

P*P,PKDATA 
R*R,RANGEI 
1*1 + 1 

* (I<,COUNT ) /LOOP 

A** DATA CONTROL VECTOR USED FOR GRAPHS IN GRAFSTAT ** 

q***************************************************** 

A ' PK= ' 

a7 4 9 4 10 pp 

a 'RANGE=' 

A 5 0 $ 4 10 pp 

C* ( ( DATA + 1 ) , ( p RANGES ) ) p RANGES 

C *, ($C) 

7 



V GUN ; I ; VAC ', BODY ; TOP ; SIDE ; VRX ; VRY ; VRZ 
q*************************************************** 
q******* SIMULATES AIR DEFENSE GUN BY CALLING ***** 
q******* APPROPRIATE SUBROUTINES ******** 

q*************************************************** 

RANCEI*\0 

1*0 

PKDATA * 10 
LOOP 3 : RANGECHEK 



34 



[9] 

CIO] 

[11] 

[ 12 ] 

[13] 

[14] 

[15] 

[16] 

[17] 

[18] 



[ 1 ] 

[ 2 ] 

[3] 

[4] 

[5] 

[ 6 ] 

[7] 

[ 8 ] 

[9] 

[ 10 ] 
[ 11 ] 
[ 12 ] 

[13] 

[14] 

[15] 

[16] 



[ 1 ] 

[ 2 ] 

[3] 
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[5! 

[7] 
[81 

[9] 

[ 10 ] 
[ 11 ] 
[ 12 ] 

[13] 

[14] 

[15] 



RADAR 

AIRCRAFT 

INTERCEPT 

MISSVECTOR 

PHVECTOR 

PROBKILL 

RANGEI+RANGEI , RTGT 
PKDATA<rPKDATA ,PK 
I>I + 1 

+(IZDATA)/LOOP3 

V 



V RANGECHEK \R\0\D\ MAXALT ; MIN ALT 

ft*************************************************** 

fi * * * * DETERMINES A/C POSITION AT TIME OF FIRING * * 

p *************************************************** 

a ' WHAT IS THE MAXIMUM ALTITUDE IN METERS? ' 

MAXALT *- 1100 

R ' WHAT IS THE MINIMUM ALTITUDE IN METERS? ' 

MINALT+90 

R+RANGE 

0+ 0 1. 5707963 UNIF 1 
+(R>MAXALT)/LOOPl 
ZF*-(MINALT ,R) UNIF 1 
+LOOP 2 

LOOP 1 : ZF+ {MIN ALT .MAXALT ) UNIF 1 
LOOP2:D+((R*2 )-(ZF*2 ))*0. 5 
XF+D*(2oO) 

YF+Dx{loO) 

V 



V RADAR ;S 2; THETAA ; R ; DELRANGE ; DTHETA ; THETAE 

p************************************************* 

p**** DETERMINES RADAR LOCATED TARGET BASED *** 
p***** ON ERRORS INDUCED BY RADAR ********* 

S2-«-0 . 00 1 2 + ( 2 GRANGE ) 

DELRANGE *- 1 NORRAND 0 6.2 
DELRANGE+DELRANGEx 1 
DTHETA+2 NORRAND ( 0 ,S2 ) 

D THE TA*r DTHETA x 1 
R*rRANGE+ DELRANGE 

THETAA+~lo(XF* ( ( (XF*2 )+ (IF* 2 ) )*0 . 5 ) ) 

thetaa+ihetaa 

THETAE *r lo (ZFiRANGE) 

THETAE*-THETAE 
YR+Rx (2o( DTHETA [ 2 ] +THETAE 
XR*-Rx ( 2 o ( DTHETA [ 2 ] +THETAE 
Z R+Rx ( 1 o ( DTHETA [ 2 ] +THETAE 

V 




2 o ( DTHETA [ 1 ] +THETAA } ) 
1 o ( DTHETA [ 1 ] +THETAA ) ) 



V AIRCRAFT \ VMIN : 4XMAX : AXMIN ; A ZMIN \AZMAX\D\0\VZ 1 
;7Z2;MIMlLr;MXALF;AXl;>lZl ;AZ2;01 

p***** A/ C PERFORMANCE DATA INPUTS ********** 

P************************************************** 

R 1 WHAT IS THE MIN ALT THE A/C IS ALLOWED TO FLY?' 
MINALT+90 

R 1 WHAT IS THE MAX ALT THE GUN WILL ENGAGE? ' 

MAXALT*- 1100 

fl ' WHAT IS THE MIN VELOCITY OF A/ C IN METERS /SEC? ' 

VMIN+150 

R ' WHAT IS THE MAX VELOCITY OF A/ C IN METERS /SEC? ' 

VMAX+ 220 

fl ' MAX A /C ACCELERATION IN Z DIRECTION IN M/SEC2 ? ' 
AZMXX+V9 

fl ' MIN A/C ACCELERATION IN Z DIRECTION IN M/SEC2 ? ' 
AZMIN*r 39 . 2 
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[16] 

[17] 

[18] 

[19] 

[ 20 ] 
[ 21 ] 
[ 22 ] 

[23] 

[24] 

[25] 

[26] 

[27] 

[28] 

[29] 

[30] 

[31] 

[32] 

[33] 

[34] 

[35] 

[36] 

[37] 

[38] 



A ' MAX ACCELERATION IN X OR Y DIRECTION IN M/SEC2 ? ' 
AXY«-0 .VxAZMAX 

q ******* ****************************************** 
fl**** DETERMINES A/ C FLIGHT PATH MANEUVER ***** 
a***** AT TIME OF FIRING *********** 

q ******************************* ****************** 

VA*- ( VMIN. VMAX ) UN IF 1 
0 <■ 0 6.2831853 UNIF 2 
7Z2«-(((7A*2)42)*0.5) 

VZ1 <r((MINALT-ZF)* (MAXALT-MINALT) )*V Z2 
VZ+CVZl t VZ2 ) UNIF 1 
D+(lVA*2)-lVZ*2 ) )*0 . 5 
VX+Id* (2O0C1] ) ) 

yy<-(cx(ioo[i] >1 

AZl+(.MINALT-ZF)*((MAXALT-MINALT)*UZMINx0.6 ) ) 

AZ2+0 .6*AZMAX 

AZ<-(AZ1,AZ2 ) UNIF 1 

01«- 0 6.2831853 £WIF 1 

AX<-AXYx(2otfl ) 

AY<-AXYx(lo01) 

7FZ<-1 NORRAND ( VZ , ( 0 . 1 , VZ ) ) 

VRX* 1 NORRAND l VX , ( 0 . 1 , VX ) ) 

VRY+l NORRAND ( VY , ( 0 . 1 , VY ) ) 



V INTERCEPT ; F12W ; BURST IME ; XF 1 ; YF1 ; ZF1 ; RPROJ ; /V ; FTN 

1] A************************************************** 

[2] a** DETERMINES PREDICTED INTERCEPT POINT BASED ** 

[3] A * * * ON LINEAR FIRE CONTROL COMPUTATIONS ****** 

[4] p*** AND THE NEWTON-RAPHSON TECHNIQUE ******* 

5] a************************************************** 
6 ] DELT+l 4 ( FI RERATE i 6 0 ) 

[7] BURST IME+BURST* (FIRERATE*6Q ) 

[8] JV4-L(BW?Sr*2) 

[9] TO+NxDELT 

[10] ZViy-e-O 

[11] LOOPl:XPI+XR+(VRXxTO )+(VRX*TAU) 

[12] ypi«-YF+(yFYxro ) + (. vry * tau ) 

[13] ZPI<rZR + (VRZxTO) + (VRZ*TAU) 

[14] RTGT+( (XPI * 2 ) + (YPI*2 ) + ( ZPI*2 ) )*0. 5 

[15] RPROJ<-(VOxTAU)*((B*TAU) + l ) 

[16] FTN<-RTGT- RPROJ 

[17] F12W«-( (XPIx7PX)+(YPIx7FY)+(ZPIx7FZ ) )*RTGT 

[18] F12W«-F1ZW- (704 (l + (BxZV4iO*2 ) ) 

[19] +(RTGT<3000)/LOOP3 

[20] AIRCRAFT 

[21] LOOP 3 :■*■(.(. \FTN)<0 .OOOD/LOOP2 

[22] TAU+TAU-(FTN*F1TN ) 

[23] +LOOP1 

[24] LOOP21+0 

V 



V MI SSV ECTOR 

I] A***************************************************** 

[2] a** CALCULATES DISTANCE FROM A/ C TO MEAN THEORETICAL * 

; 3 ] a** INTERCEPT POINT FOR EACH PROJECTILE ******* 

4] a***************************************************** 

[5] APN+x 0 

[6] 7IMPiV«-t0 

[7] CPAJV-M0 

[8] N + 1 

[9] TAUQ+ 0.1 

[10] LOOP : MISSDIST 

II] A***************************************************** 
112] a** COLLECTS STATISTICS ON BIS , PRESENTED AREA, *** 

[13] a** AND BURST VELOCITY . **** 

[14] a** USED TO DETERMIN PH IT AND PRILL ******* 

15] a***************************************************** 
[16] CPAN+CPAN ,CPA 
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[17] 

[18] 

[19] 

[ 20 ] 
[ 21 " 
[ 22 ] 



APDATA 
APN+APN ,AP 
VIMPN+VIMPN , V IMP 
N+N+l 

(NZBURST ) /LOOP 
-*■0 
V 



:i] 

! 2 ] 

!3] 

’ 4 " 
!5] 
[ 6 ] 

[7] 

[ 8 ] 

[9] 

[ 10 ] 
[ 11 ] 
:i2] 
:i3] 
:i4] 

[15] 

[16] 

[17] 

:is] 

[19] 

[ 20 ] 
. 21 ] 
! 22 ] 
;23] 
.24] 
!25] 
!26] 
.27] 
!28] 
[29] 
’.3 0] 
!3 1] 
!32] 
[33] 
]34] 
!35] 
!36] 
!37] 
!38] 
!3 9] 



7 MISSDIST ; A: FI ; F 2 ; F3 ; XA ; YA ; ZA ;F1TN; F2TN ; BT 2 ; BT 
; XAP ;YAP; ZAI?;TQ 

P************************************************ 

p **** DETERMINES TIME AT WHICH CLOSEST POINT ** 
p * * * * * OF APPROACH OCCURS ****** 

q************************************************ 

VPX+(XPI*RTGT)*V 0 

vpy+IypurtgtUvo 

VPZ<r(ZPIiRTCT)xVO 
LOOPl:BT+(B*TAU)+l 
BT2<rBT*2 
XP+(VPXxTAU)*BT 
YP+lVPYxTAU}*BT 
ZP+lVPZxTAUUBT 
T0<-(N-1 )*DELT 
XA+XF+ (VX*T0 )+( 

+ (AX (TAU*2) 2 
YA+YF+(VY*T0 )+( 

+ (AY (TAU*2) 2 
ZA+ZF+(VZ*T0)+( 

+ (AZ (TAU*2) 2 
XAP+XA-XP 
YAP+YA-YP 
ZAP+ZA-ZP 

Fll*-XAPx (VX+ (AXxTO ) + (AX*TAU 
F12<rYAPx(VY+(AYxT0 ) + (AYxTAU 
F13+ZAP*(VZ+(AZxT0 )+(az*tau 
F14<-( (XAP*2 ) + (YAP*2 )+(ZAP*2 
F 1 TN*- (F11+F12+F13 ) 

A<-BT* 3 

F21<-( (F11*XAP)*2 
F22<-( (F12*YAP)*2 
F23<-( (FI 3 * ZAP )*2 
F2TN+F21+F22+F23 
■* ( (TAU-TAUO )<0 .0001 )/LOOP2 
TAUO+TAU 

TAU+TAU- (F1TN+F2TN ) 
fi J-e-J + 1 

fi-*( (LOOP2*\I)*30 )/LOOP2 
+LOOP 1 

LOOP2:CPA+Fin 

TAUQ+Q 

VP+VOi ( ( (B*TAU)+ 1 )*2 ) 

RPROJ<-(VOxTAU)*( (B*TAU)+1 ) 

V 



AX*2 )x (T 0*2 ) ) + ( (VX+(.AXxT0 ) )*TAU) 
AY *2 )x(ro*2 )) + ( (VY+(.AYxT0 ) )*TAU) 
AZ*2 )*(T0*2) )+((VZ+(AZ*T0))xTAU) 



\-(VPX*BT 2 
)-(VPY*BT 2 
\-(VPZ*BT 2 
)* 0. 5 




2xVPXxB*A 
2 xVPY*B*A 
2*VPZ*B*A 



\*xap: 

lx YAP 
IxZAP. 



V APDATA ; VXIMP ; VY IMP ; VZIMP ; API ; AP2 

1] p************************************************* 

1 2 ] fl ** COMPUTES PRESENTED AREA OF THE AIRCRAFT ** 

3] q ******************** ***************************** 
’.4] VXIMP+VX-(VP*XPiRPROJ) 

:5] VYIMP+VY-lVPxYP*RPROJ} 

!6] VZIMP+VZ- ( VP*ZPiRPROJ ) 

[ 7 ] VIMP+ ( ( VXIMP* 2 ) + ( VYIMP* 2 ) + ( VZIMP* 2 ) ) *0 . 5 
] 8 ] + (VIMP<>VP)/LOOP , ,,, 

[9] APl^ASIDE- ((VIMP-VP)x (ASIDE- AFRT)*(\ (7A))> 

[10] AP2+ATOP-((VIMP-VP)x(ATOP-AFRT)*( | (VA ) ) ) 

"11] +LOOP2 

: 1 2 ] LOOP : AP1 + ( ( V IMP -VP ) x (ASIDE-AFRT )* ((I (VA ) ) ) ) +ASJPF 
AP2+ ( (VIMP-VP)* (ATOP-AFRT )* ( (\ (VA))))+ATOP 



:i3 

:i4 



LOOP2:AP+(APl,AP2) UN IF 1 
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V 




V PHVECTOR 

q************************************************ ★* ★ 
fl** STORES PROBABILITY OF KILL FOR EACH ROUND **** 

q**rt^*****rt^rt^***rt****rtrt****rt***rt*******rt*rt********* 

PH IT * t 0 
N*pCPAN 
I* 1 

L00P:CPA*CPANZI1 
PROBHIT 
PHIT*PHIT ,PH 
1*1 + 1 

* (I£N ) /LOOP 

V 



V PROBHIT ; I ;N 

[1] A************************************************* 

[2] o*** CALCULATES PROB OF HIT FOR A SINGLE ROUND ** 

[3] B************************************************* 

[4] R*RANGE 

[5] A«-((6.2831x(fl*2 )*(DISP* 2) )+AP) 

[6] B*\ -3 . 14 21* ( CP A* 2 ) ) 

[7] C<-(2x(R*2')x{DISP*2))+AP 

[8] D*B+C 

[ 9 ] PH+{AP*A)*i*D) 

V 



Cl] 

C 2 ] 
C3] 
C4] 
[5] 

C 6 ] 
[7] 

cs; 
cs: 
cio] 
!iij 
C 1 2 ] 
C 1 3 ] 
C 1 4 ] 
Cl 5] 
C 1 6 ] 
C 1 7 ] 
C18] 



V PROBKILL ;I 

q**************************************************** 
o * * DETERMINES PROBABILITY FOR A GIVEN BURST * * * * 

B**************************************************** 

Y* 0 0 0 1 
1*1 

L00P1 : VIMP*VIMPNZn 

a** MINIMUM IMPACT VELOCITY THAT CAUSES DAMAGE **** 



B**************************************************** 

* ( VIMP<, 500 ) /STOP 
AVDATA 
PH*PHITII 1 
*(PH<1E Q)/L00P2 
*( (+/X)=0 )/L00P 2 



MARKOV 

L00P2:I*I+1 



*(IZBURST)/L00P 1 
ST0PiPK*(,Y) C13 
7 



V A ; VMIN \MAXV \AVMAX\A VMIN 

Cl] B***************************************************** 

C 2 ] n** BUILDS TRANSITION MATRIX FOR VELOCITY AT BURST *** 

C 3 ] a***************************************************** 

C 4 ] 0 'AV DATA FOR MIN VELOCITY EXPRESSED AS A VECTOR ? ' 

C 5] o AVMIN *1 0.11030.11030.082300.889700.02800 
0.8897 0.028 0 0 0 0.8617 

C6] a AVMIN *1 0.1408 0.1408 0.0988 0 0.8592 0 0.0421 0 0 
0.8592 0.0421 0 0 0 0.817 

C7] AVMIN * 10.20180.20180.131700.798200.070100 
0.7982 0.0701 0 0 0 0.7281 

C8] fl AVMIN* 1 0.26 28 0.26 28 0.1646 00.73 72 00.098200 
0.7372 0.0982 0 0 0 0.639 

C9] n 'AV DATA FOR MAX VELOCITY EXPRESSED ASA VECTOR ? ' 

CIO] aAVMAX*l 0.2073 0.2073 0.1372 0 0.7297 0 0.0701 0 0 
0.7927 0.0701 0 0 0 0.7226 

Cll] aAVMAX* 1 0.23 41 0.23410.1500.76 5900.0 84100 
0.7659 0.0841 0 0 0 0.6818 

C 12] AVMAX* 1 0.2018 0.2018 0.1728 0 0. 7982 0 0. 1122 0 
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[13] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[ 20 ] 
[ 21 ] 
[ 22 ] 



0 . 7982 0 . 1122 0 0 0 0.6028 

qAVMAX a 1 0.3165 0.3165 0.1902 0 0.6835 0 0.1262 0 0 
0.6835 0.1262 0 0 0 0.5573 
fi ' WHAT IS THE MINIMUM VELOCITY? ' 

VMIN+50 

+ (VBUST<, 100 ) /STOP 

fi ' WHAT IS THE MAXIMUM VELOCITY? ' 

MAXV+VQ+VMAX 

Xa ( ( < VBUST-VMIN ) x (AVMAX-AVMIN ) ) * ( MAX V - VMIN ) ) +A VMIN 
+GO 

STOP : X+0 
G0:+0 
V 



V MARKOV: DIMENSION: AVMATRIX ; IDMATRIX ; TRANSMATRIX 
; YNEW ; COUNT ;PM;L ;N;X1 

[1] A**************************************************** 

[2] A**** BUILDS LETHALITY MATRIX AND TRANSITIONS ****** 

[3] a**** IT TO GIVE PROBABILITY OF KILL ************** 

[4] a**************************************************** 

[5] a ' INPUT VULNERABLE AREA DATA AS AN ARRAY ' 

[6] Xl+X 

[7] La pXl 

[ 8 ] DIMENSIONAL * 0 . 5 

[ 9 ] A VMATRIX+ (DIMENSION , DIMENSION ) pXl 

[10] A ' INPUT THE INITIAL STATE VECTOR r 

[11] Al*0 0 0 1 

[12] a ' WHAT IS THE PROBABILITY OF HIT ' 

[13] fi PH+ti 

[14] A VMA TRIX+A VMA TRIX x PH 

[15] PM+l-PH 

[ 16 ] IDMATRIXa ( ( i DIMENSION ) « . x ( i DIMENSION ) ) 

[17] IDMATRIX a ( ( ( DIMENSION , DIMENSION ) p ( ( i DIMENSION ) * 2 ) ) 

= IDMATRIX) 

[18] TRANSMATRIXa( (PM*IDMATRIX)+AVMATRIX) 

[19] A ' THE NUMBER OF MATRIX TRANSITION ' 

[20] Na l 

[21] GO:Ya( DIMENSION, DIMENSION ) p Y 

[22] YNEWaTRANSMATRIXxY 

[23] YNEWa+/(§YNEW) 

[24] YNEWa(i,DIMENSION)pYNEW 

[25] Ya^YNEW 

[26] COUNT aN- 1 

[27] NaCOUNT 

[28] a(cOUNT>0)/GO 

V 




V CUMKILL 

Q'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

A** GIVES CUMULATIVE PROBABILITY OF KILL FOR **** 
a** A TARGET THAT HAS BEEN ENGAGED WITH MORE **** 

A*** THAN ONE BURST ******** 

q************************************************* 

A$(((DATA+l),(pRANGES))pP) 

Al-C 

S<3) 




[ 1 ] 

[ 2 ] 

[3] 

[4] 

[5] 

[ 6 ] 
[7] 



V RaAB UNIF N\A\B 

Q'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

A* GENERATES N UNIFORM RANDOM NUMBERS ********** 

A*** ON INTERVAL (A ,B ) ********************** 
a******** * * *************************************** 

aaabl i] 

fl<-4fl[2] 

+ )x (*2147483647 )x ?tfp 2 14 748 3 646 
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V 



V R+N UNIRAND B 

Cl] fi************************************************* 

[ 2 ] p * * GENERATES N PSEUDO UNIFORM RANDOM NUMBERS * * * 

[3 ] q************************************************* 

[4] £f(B[l]-B[2] ) + ( OV?100000000 )x2x£[2] )*100000000 

V 

V Z+/V NORRAND P;S;I;T 

"1] q********************************************* 

12] q*** GENERATES N NORMAL RANDOM NUMBERS ***** 

'3] q********************************************* 

|4] P«-2+P,l 

[5] Z«-(lV)pO 

[6] It - 1 

[7] F10:T+2 UNIRAND 0.5 0.5 

[8] ZV(2xr)-l 

[9] 5^(r[l]*2)+r[2]*2 
110] -*F10xiS>l 

[11] Z[J]^P[l]+P[2]x(Z’[l]x(( 2x®S)*S)*0.5) 

[12] ->-F10xx (J<-I + l )<,N 

V 



CORRECTIONS MADE TO KEELING'S PROGRAM IN : 
Subroutine AIRCRAFT (36, 37, 3S) 

Subroutine MISSVECTOR (9) 

Subroutine MISSDIST (5, 30, 31, 37) 

Subroutine PROBHIT (6) 
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APPENDIX B 

APL COMPUTER PROGRAM 



1. NEW PROGRAM 



1 

$ 
£13 

7 

8 
[ 9 ] 
[10 
'll 
[ 12 ] 

13 ] 

141 

15 ' 

16 

17 

[18 

[19 

[ 20 ] 

[21* 

[22 

[23 

'24 

r Pi 

L2o] 

27 ] 

28 ] 
[ 29 ] 
' 30 " 

31 

32 

[ 33 ] 

' 34 ' 

35 ] 

36 ' 

[37 

[ 38 ] 

[ 39 ] 
[40 
[41 
[42 

[ 43 ] 

[ 44 ] 

[ 45 ] 

[ 46 ] 
[47 ~ 
[ 48 ] 



V DRIVER', COUNT', I 

A*rt*********rt*****rt*******rtrt**rt*rt******rt*******rt**rt 

o******* DATA INPUT ******************** 

p************************************************** 

a « HOW MANY DATA POINTS DO YOU WANT PER RANGE?' 

DATA * 25 
DAT A* DATA - 1 

a 'WHAT IS THE GUN FIRING RATE IN ROUNDS PER MINUTE?' 
FIRERATE+6 0 0 

fi 'WHAT IS THE MUZZLE VELOCITY?' 

70+1005 

A 'WHAT IS THE BURST SIZE?' 

BURST* 5 

A 'WHAT IS THE ANGULAR DISPERSION IN RADIANS ? 1 
PISP+0.001 

A * WHAT IS THE DRAG FACTOR B=CD*V 0*0.5?' 
fl+0.163 

a 'WHAT ARE THE RANGES YOU WANT TO LOOK AT?' 

RANGES * 1900 1300 700 100 
a 'WHAT ARE THE AREAS OF THE AIRCRAFT? ' 

AFRT * 5 
ASIDE * 30 
ATOP * 50 
RL* 5 

PK/H* 0 .25 

A*************************************************** 

p*** PROGRAM EXECUTION AND DATA COLLECTION ******* 
p*** P IS PROBABILITY OF KILL ****** 

a*** R SLANT IS RANGE OF TARGET AT INTERCEPT *** 

p *************************************************** 

R *\ 0 
P+iO 
1*1 

COUNT* oRANGES 

LOOP : RANGE*RANGES [ J] 

GUN 

P*P ,PKDATA 
R*R , RANGE I 
1 * 1+1 

*(I£COUNT) /LOOP 

A***************************************************** 

a** DATA CONTROL VECTOR USED FOR GRAPHS IN GRAF ST AT ** 

p ***************************************************** 

a 'PK=' 

a 7 4 ? 2 10 p P 
a ' RANGE- ' 

A 5 0 5 4 10 p R 

nC*C(DATA + l ) , (pRANGES) )p RANGES 
aC+ , (<5C ) 

7 



V GUN ; I •, VAC ', BODY ', TOP ', SIDE ', VRX ; VRY ; VRZ 

A*************************************************** 

p******* SIMULATES AIR DEFENSE GUN BY CALLING ***** 
p******* APPROPRIATE SUBROUTINES ******** 

A ********* ****************************************** 

RANGE I * \ 0 
1+0 
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'.71 

:s] 

;g] 

:io] 

:ii] 

:i2] 

:i3] 

:i4] 

:i5] 

:i6] 

:i7] 

:i8] 



PKDATA +\ 0 

L00P3 :RANCECHEK 

RADAR 

AIRCRAFT 

INTERCEPT 

MISSVECTOR 

PKVECTOR 

BURSTKILL 

RANGE I GRANGE I , RTGT 
PKDATA+PKDATA ,PK 
I+I + l 

+(IZDATA)/L00P3 

V 



V RANGECHEK ;R;0 ; D ; MAX ALT ; MIN ALT 

1] A*************************************************** 

[2] a**** DETERMINES A/C POSITION AT TIME OF FIRING ** 

[4] a 'WHAT IS THE MAXIMUM ALTITUDE IN METERS?' 

[5] MAXALT + 1100 

!6] a 'WHAT IS THE MINIMUM ALTITUDE IN METERS?' 

!7] MI/W1LZV9 0 
!8 3 P+P^l/VCF 

!9] ^ 0 1. 5707963 UNIF 1 

:i0] +(R>MAXALT)/LOOP 1 
:il] ZF*-(MINALT ,R) UNIF 1 
:i23 +LOOP2 

:i3] L00Pl:ZF+(AmV4£Z\AMX40r) UNIF 1 
:i4] P00P2 :£«-( 0?*2 )- (ZF*2 ) )*0 . 5 
:i5] XF*-D* ( 2 oO{ 

:i6] YF<rD*lloO) 

V 



V RADAR ;S2; THETAA ; R ; DELRANGE ; ; THETAE 

13 A************************************************* 

[23 a**** DETERMINES RADAR LOCATED TARGET BASED *** 

33 a * * * * * ON ERRORS INDUCED BY RADAR ********* 

"43 A************************************************* 
153 S2-<-0.0012 + (2+RANGE ) 

6 3 PFLPd/VCF-^l NORRAND 0 6.2 

7 3 PFLF4/V0F+FFLP4W0FX 0 . 1 
[83 DTHETA+2 NORRANDiO ,S2 ) 

[9 3 DrPETd-f-PrffETyixO • 1 
[103 R+RANGE+DELRANGE 

[113 THETAA+ lo(XF*(( (XF* 2 )+(YF* 2 ))*0. 5 )) 

[123 

[13 3 THETAE+ lo(ZF*F4M7F) 

[143 THETAE+THETAE 

[153 YR+R x(2o (DTHETA [ 2 3 +PP ET4E ) ) x ( 2 o (DFPFZ71 [ 1 3 +THETAA ) ) 
[16 3 XP+Px ( 2 o ( DTHETA [ 2 3 +THETAE } ) x ( l o {DTHETA [ 1 ] +THETAA ) ) 
[173 ZR+Rx{lo{DTHETAZ21+THETAE)) 

V 



V AJPCPAFF : VMIN : AXMAX : AXMIN ; 4 ZMIrt ; A ZMAX ; D ; 0 ; 7Z1 
;7Z2 ;MIMLr;MAXA£r;AXY;AZl ;4Z2 ;01 
.13 A************************************************** 
! 2 1 a***** 4 /C PERFORMANCE DATA INPUTS ********** 

3 J fl************************************************** 
[43 a ' J/PAr IS TPFMIJV ALT rffFd/C IS ALLOWED TO FLY?' 

[53 MINALT+ go 

[6 3 a ' WHAT IS THE MAX ALT THE GUN WILL ENGAGE? ' 

[73 MAXALT+1100 

[8 3 a 'WHAT IS THE MIN VELOCITY OF A/C IN METERS/ SEC? ' 

[9 3 7MIJV+150 

[10 3 a 'WHAT IS THE MAX VELOCITY OF A/C IN METERS / SEC? ' 

[113 VMAX+ 220 

[12 3 a ' MAX A/C ACCELERATION IN Z DIRECTION IN M/SEC2? ' 

[133 AZMAX+ 49 
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[14] 

[15] 

[16] 
[17" 
[18] 
[19] 
[ 20 " 
[ 21 ] 
[22 
[23: 
[24 
[25 
[26: 

[27] 

[28] 
[29" 
[30] 
[31" 
[32 
[33: 
[34 
[ 35 : 
[36 
[37 
[38 



a ' MIN A/C ACCELERATION IN Z DIRECTION IN M/SEC2 ? ' 
AZMIN+3 g . 2 

a ' MAX ACCELERATION IN XORY DIRECTION IN M/SEC2 ? ' 
AXY+O .UxAZMAX 

P ************************************************* 

a**** DETERMINES A/C FLIGHT PATH MANEUVER ***** 
p***** AT TIME OF FIRING *********** 

P ************************************************* 

VA+ ( VMAX, VMAX ) UN IF 1 
0*- 06.2831853 UN IF 2 
VZ2+(((VA*2)*2)*Q .5) 

VZl + ( (MINALT-ZF)i (MAXALT-MINALT) )*VZ2 
VZ+(VZ2.VZ2 ) UNIF 1 
D+((VA*2')-IVZ*2))*0.S 
VX+lD*(2oOlll )) 

VY+(D* (lo0[l] ) ) 

AZ1+(MINALT-ZF)* ( (MAXALT-MINALT)* (AZMINxQ . 6 ) ) 
AZ2+0.6xAZMAX 
AZ«-(AZ2,AZ2) CWIF 1 
01+ 0 6.2831853 CWIF 1 
AX+AXYx (2o01 
AY+AXYx(lo01 
VRZ+1 NORRAND 
VRX+1 NORRAND 
VRY+1 NORRAND 
V 




[ 1 ] 

[ 2 ] 

[3] 

'4' 

[5] 

[ 6 ] 

[7] 

[ 8 ] 

[9] 

[ 10 ] 
[ 11 ] 
[ 12 ] 

[13] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[ 20 ] 
[ 21 ] 
[ 22 ] 

[23] 

[24] 



V INTERCEPT ; F1ZW ; BURST IME ; XF1 ; YF1 ; ZF1 ; RPROJ ; N ; FTN 



P ************************************************** 

A * * DETERMINES PREDICTED INTERCEPT POINT BASED * * 
a * * * ON LINEAR FIRE CONTROL COMPUTATIONS ****** 
a*** AND THE NEWTON -RAPHSON TECHNIQUE ******* 

P************************************************** 



DELT+1* (FIRERATE*60 ) 

BURSTIME+BURST* (FIRERATE*60 ) 

N+L (BURST *2 ) 

TO+NxDELT 
TAU+ 0 

LOO PI : XPI+XR* ( VRX*T0 )+(VRX*TAU) 
YPI+YR+(VRY*TQ )+(VRY*TAU) 

ZPI+ZR+lVRZxTQ )*(VRZxTAU) 

RTCT+ ( (XPI*2 )+(YPI* 2 )*(ZPI*2 ) )* 0. 5 
RPROJ+(VQ*TAU)*( (B*TAU)+1 ) 

FTN<-RTGT- RPROJ 

FlTN+( (XPI*VRX)+ (YPI*VRY )+(ZPI*VRZ ) )*RTGT 
F1TN+F1TN- ( VO * ( 1 + ( B *TA U ) * 2 ) ) 

+(RTGT< 3000 )/LOOP 3 
AIRCRAFT 

LOOP 3: + (( \FTN)<0 .OOOD/LOOP2 
TAU+TAU- (FTN*F1TN ) 

+LOOP 1 
LOOP2X+Q 
V 



V MISSVECTOR 

: 2 ] a** CALCULATES DISTANCE FROM A/C TO MEAN THEORETICAL * 

!3] a** INTERCEPT POINT FOR EACH PROJECTILE ******* 

4] R ********************** ******************************* 

! 5] APN+x 0 

!6] VIMPN+x 0 

[ 7 ] CPAN + \ 0 

! 8 ] N+l 

!9] TAU 0+0.1 

[10] LOOP : MISSDIST 

11] p***************************************************** 

[12] p** COLLECTS STATISTICS ON BIS , PRESENTED AREA , *** 

;i3] a** AND BURST VELOCITY . **** 

'14] a** USED TO DETERMINE PKSS AND PR/ BURST ******* 
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[15] A***************************************************** 

[16] CPAN+CPAN , CPA 

[17] APDATA 

[18] APN+APN.AP 

[19] VIMPN+VIMPN , VIMP 

[20] N*-N + 1 

[21] +(N£BURST ) /LOOP 

[22] ->-0 

V 



:i] 

E 13 

Z4] 

[5] 

7 

[ 8 ] 

!9] 

;i°] 

.ill 

:i2] 

:i3] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[ 20 ] 
[ 21 ] 
[ 22 ] 

[23] 

[24] 

[25] 

[26] 

[27] 

[28] 

[29] 

[30] 

[31] 

[32] 

[33] 

[34] 

[35] 

[36] 

[37] 

[38] 

[39] 



V MISSDIST: A: FI ; F2 ; F 3 ; XA ; YA ; ZA ; F1TN ; F2TN ; BT 2 ; BT 
;XAP’,YAP;ZAP;TQ 

A************************************************ 

a * * * * DETERMINES TIME A T WHICH CLOSEST POINT * * 
a * * * * * OF APPROACH OCCURS ****** 

a ************************************************ 
VPX+ ( XPURTGT ) x 70 
VPY<r(YPI*RTCT}xVO 
VPZ+(ZPIiRTGT)xV 0 
LOOP1 :BT<-(B*TAU)+1 
BT2*-BT*2 
XP*- ( VPXxTAU ) *BT 

yp+IvpyxtauUbt 

ZP+- ( VPZxTAU )*BT 
TO+lN-l )*DELT 
XA+XF+(VX*T 0 )+( 

+ (AX XTAU*2) 2 
YA<rYF+(VY*T 0 ) + ( 

+ (AY XTAU* 2 
ZA+ZF+(VZxT 0 
+ (AZ XlAU* 2 
XAP+XA -XP 
YAP+YA-YP 
ZAP+ZA-ZP 
Fll+XAPxiVX* 

F12+YAPx(VY+ 

F13+ZAP* XVZ* 



AX* 2 )x(ro*2 ) ) + ( (7X+UXXY0) )*TAU) 
AY *2 )x(Y0*2 ) )+( (VY+(AY*TO ) )*TAU ) 
AZ*2 )*(T0*2 ) )+( (VZ+(AZxTO))*TAU) 



AX*TAU] 
AY*TAU , 
AZxIAU . 




AX*TO )+ 

AY*T 0 )+ 

_ .AZxro )+ 

F14«-( (XAP*2 )+(YAP*2) + (ZAP*2 
F1TN*-(f11 +F12 +F13 ) 

A<rBT*3 

F21<rUFll*XAP)*2 
F 22-K (F12*YAP)*2 
F23-K (F13 *ZAP )*2 
F2TN+F2 1 +F2 2 +F2 3 
->■( (TAU-TAUO )<0 .0001 )/LOOP 2 
TAUO+TAU 

TAU+TAU- (F1TN*F2TN') 

A I-s-I + 1 

a-»( (,LOOP2*\I)Z30)/LOOP2 
+LOOP 1 

L00P2:CPA+FV* 

TAU 0+0 

VP+V 0*( ( ( B*TAU)+1 )*2 ) 
RPROJ+XVQxTAU)*((BxTAU)+l ) 

V 



-(7PX+Sr2 
-\VPY*BT2 
-lVPZ*BT2 
)*0 . 5 



2 xVPX*B*A, 
2 *VPY*B*A 
2 *VPZ*B*A 



>*xap: 

)xYAP 

)*ZAP. 



V APDATA ; VXIMP ; VYIMP ; VZIMP ; API ; AP2 

[I] A************************************************* 

[ 2 ] a * * COMPUTES PRESENTED AREA OF THE AIRCRAFT * * 

[3] R************************************************* 

[4] VXIMP+VX-(VPxXP*RPROJ) 

[5] VYIMP<-VY-(.VPxYP*RPROJ ) 

[6] VZIMP<rVZ- {VP*ZP*RPROJ ) 

[ 7 ] VIMP*- ( ( VXIMP* 2 ) + ( VYIMP* 2 ) + ( VZIMP* 2 ) ) *0 . 5 

[8] ->-(7IMPs7P)/L00P 

[9] APl^ASiTE-((7IMP-7P)x(ASIPF-AFPr)+(l (7A))) 

[10] AP2«-ArOP- ( (7IMP-7P)x (ArtfP-AFFP)* ( | (7A ) ) ) 

[II] +L00P2 

[12] LOOP : API-*- ( ( 7IMP- 7P ) x ( AS IDE-AFRT ) * ( ( | ( VA ) ) ) ) +ASJDF 



44 



[13] AP2<-( (VIMP-VP)x (ATOP-AFRT)* ( ( | ( VA ) ) ) )+ATOP 

[14] LOOP 2 :AP+(.AP 1 ,AP 2 ) UNIF 1 



V PKyPCTtfF; I;N 

1] fl**************************************************** 

12] n** STORES PROBABILITY OF KILL FOR EACH ROUND ** 

3] a**************************************************** 

[4] PKSSN+\ 0 

[5] N+pCPAN 

[6] 1+1 

[7] LOOP:CPA+CPANin 

8] A**************************************************** 

[9] A*** OPTIONAL SELECTION FOR TYPE OF FUZE ********** 

[10] A * * * FOR PROXIMITY FUZE. USE PROXFUZE ********* 

111] a*** FOR CONTACT FUZE , USE CONTFUZE ************** 

12] a**************************************************** 

[13] PROXFUZE 

[14] ftCONTFUZE 

[15] PKSSN+PKSSN ,PKSS 

[16] I+J+l 

[17] ■* (I<,N ) / LOOP 

V 




V PROXFUZE ;A;B-,C;DiR;RL',R 0 ;AL;PF 

A************************************************** 

a** PROXIMITY FUZED WARHEAD STORES PROBABILITY ** 

A*** OF KILL FOR SINGLE SHOT ********** 

A************************************************** 



[ 5 ] R+RANGE 

[6] PF+ 1 

[7] RL*- 5 

18] RQ*-( 1 . 2*RL ) 

[9] 4+(P0*2)xPF 

[10] B*- ( 2 x (DISP*2 )x(F*2)) + (F0*2) 

[ 11 ] 04*5 

[12] ZM 1 x (CPA *2 ))•*■(( 2x (DISP*2 )x(F*2))+(P0*2)) 

[13] PFSS+Cx(*P) 

V 



V CONTFUZE;A;BiC;DiR 

[1] A*********************************************** 

[2] a * * * CONTACT FUZED WARHEAD CALCULATES ******** 

[3] ft*** PROBABILITY OF KILL FOR SINGLE SHOT ***** 

[4] a**** AV-PK/H^AP, ASSUMED PK/H-0 . 25 ********** 

[5] A*********************************************** 

[6] R+RANGE 

[7] 4<-0.25x4P 

[8] B+((2x3.142)x(PISP*2)x(P*2)) + UP*2) 

[9] C+A*B 

[10] P+( - 3 . 142x (CPA* 2 ) )*S 

[11] PKSS+C*(*D) 

V 



V BURSTKILL;AK;BK 

[1] q************************************************ 

[2] ft*** PROBABILITY OF KILL GIVEN BURST *********** 

[3] R*********************************************** 

[4] ftAK+ ( (pRANGES ) , (DATA + lDpPKSS 

[5] AK+PKSSN 

[6] BK+l-AK 

[7] CK+x/BK 

[8] EK+l-CK 

[9] PK<rEK 

[ 10 ] R************************************************ 

[11] a**** PROBABILITY OF KILL ENCOUNTER ********** 

[12] a****** ASSUMED PD*PF=1 ******************* 

[13] a***** PK/E EQUAL TO PKBURST ******************* 
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[14] A ************************************************ 
V 



V R+AB UNIFN\A\B 

1] a************************************************* 

12] A * GENERATES N UNIFORM RANDOM NUMBERS ********** 

13] a*** ON INTERVAL (.A ,B ) ********************** 

"4] A************************************************* 

[5] A+AB [1] 

[6] B+ABl 2] 

[7] + )x ( + 214748 3647 ) x ?tfp 2 14748 3 646 

V 



V R+N UNIRAND B 

1] Q************************************************* 

12] A * * GENERATES N PSEUDO UNIFORM RANDOM NUMBERS *** 

3] a************************************************* 
[4] B-«-(B[l] -B[2] ) + ((/V?100000000)x2xf}[2] )*100000000 

V 




V Z+N NORRAND P;S;I;T 

Q********************************************* 

a*** GENERATES N NORMAL RANDOM NUMBERS ***** 

a********************************************* 

P-e-2 tP.l 
Z«-(AOpO 
I *- 1 

FI 0 : T+2 UNIRAND 0.50.5 
F-e- ( 2 *T ) - 1 
S<Kr[l]*2 )+Z’[2 ] *2 
^•FlOxxSSl 

Z[J]«-P[l]+PC2]x(r[l]x( ( 2x®S)*S)*0. 5 ) 

-►FlOx i (J-«-I+l )£/V 

V 
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APPENDIX C 
SYMBOL NOTATION 



SYMBOL 


DEFINITION 


UNITS 


jjp jJ; ;J; .J. .1* ;Jj ><- A A ;J 


‘r. :J? & sjt * j}j >jt # # # rjr >}c jjs $ >’/. t'f. jfc * # # jfs * jjs * >Js * * >Js # jfc * £ 


A- ^ A A A -A* A 


,A Trt 


area of aircraft front 


2 

meters 


Atop 


area of aircraft top 


i 

meters" 


Aside 


area of aircraft side 


meters" 


A max 


max aircraft acceleration 


•> 

meters sec" 


A • 

4 nun 


min aircraft acceleration 


meters, sec" 


Ap 


aircraft presented area 


2 

meters" 


A\ 


aircraft vulnerable area 


-> 

meters" 


Azl’Az2 


range of A z 


■> 

meters sec" 


A v A y ,A z 


acceleration of aircraft 


meters/sec" 


C d 


coefficient of drag 


1 meters 


g 


gravitational constant 


meters sec" 


GD 


gravity drop 


meters 


LM 


lethality matrix 


n/a 


P H 


probability of aircraft hit 


n/a 


P h 


probability of component hit 


n, a 


Pk 


probability of kill 


n/a 


P ke 


probability of kill engagement 


n, a 


P k|h 


probability of kill given hit 


n/a 


R 


slant range 


meters 


R P 


projectile range 


meters 


t 


time 


sec 


V a 


aircraft total velocity 


meters ’ sec 


^ bust 


burst velocity 


meters, sec 


V 0 


muzzle velocity 


meters sec 
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V P 


projectile velocity 


meters sec 


vw 


projectile velocity component 


meters sec 


V V V 
rx’ ry» rz 


radar aircraft velocity 


meters sec 


V \ ,V y ,V z 


aircraft velocity components 


meters sec 


V zl- V z2 


range of V 


meters sec 


w i 


uncorrected wind 


meters sec 


V Y 7 

1 a 1 ci 


aircraft position at time t 


meters 


Xf-Yf-Zf 


aircraft position at firing 


meters 


x y z 
'V p" p 


projectile location 


meters 


X V z 

‘V r r 


radar located target 


meters 


a dx ,a dy 


angular ballistic dispersion 


radians 


AX p 


projectile displacement 


meters 


C R 


range error 


meters 


e 0e- c 0a 


azimuth and elevation error 


radian 


Mbx’^by 


linear mean bias 


meters 


(co e .co a) 


gun azimuth and elevation 


milliradians 


P 


projectile linear range 


meters 


^dx’^dv 


linear ballistic dispersion 


meters 


t 


projectile flight time 


sec 


0 bv 0 b y 


angular mean bias 


radians 


0 e .0 a 


target azimuth and elevation 


milliradians 
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APPENDIX D 
EQUATIONS AND DATA 



I. GOV ERNING EQUATIONS 

a. The Siacci Formula 

This equation used here to determine the velocity of the projectile at a point in time is 
v p = V 0 (C d v o 1;2 T+l) 2 (eqnD.l) 

Where these two equation can be seen [Ref. 1]. 

• Cj is the coefficient of drag of the projectile. 

• Vp is the velocity of the projectile. 

• Vq is the muzzle velocity. 

• x is the time of flight of the projectile. 

b. Range of the Projectile 

The formula for determine the range of the projectile at a point in time is taken as 
R g = V 0 T; ( C d V 0 I,2r+1 ) (eqnD.2) 



The derivation of these formula can be found in [Ref. 4] 

c. Gravity Drop 

GD = (-1, 6)g T 2 (l + 2(V p /V 0 ) l/2 ) (eqn D.3) 



2. DATA 

This study looked at 40 mm gun system, and firing high explosive round with 
both a contact fuze and a proximity fuze. The unclassified system characteristics were 
obtained from [Ref. 5] 
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All references to vulnerable area data for specific air threats to these guns are classified, 
therefore aircraft probability of kill given a hit and warhead lethal radius is entirely 
synthetic. Vulnerable area data for many different aircrafts and many different gun 
systems is readily available for personnel with appropriate clearances. The optimal 
angular ballistic dispersion was determined using the computer simulation and the 
coefficient of drag for the projectiles are assumed, the other data was from Jane's 
Weapon System. The aircraft's vulnerability measure rj is assumed as discussed in 
Chapt 2, and Ay is determined by using the probability the aircraft is killed given the 
random hit on the aircraft P K H as discussed in Chapt 3. 



TABLE 1 

40 MM GUN CHARACTERISTICS 


parameters 


unit 


data 


firing rate 


rds/min 


600 


muzzle velocity 


m/sec 


1005 


coefficient of drag 


1/sec 


0. 00515 


maximum range 


km 


12 


angular dispersion 


mrad 


1 . 0 
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APPENDIX E 

EMPIRICAL COMPARISON OF MARGINAL DISTRIBUTIONS 

1. TYPE OF FUZE 

a. Proximity Fuzed Warhead 



TABLE 2 

PROXIMITY FUZED WARHEAD 



PARAMETER VALUES 
POPULATIONS 
WEIGHTS 
X LABEL 
Y LABEL 



(100); (700); (1300); (1900) 
RG14 ;RG13 ;RG12 ;RG11 
0 

RANGE IN METERS 
PROBABILITY OF KILL 



POPULATION NO. OF PERCENTILES 



NUMBER 


POINTS 


Y MEAN 


Y STD DEV 


0.25 


0.5 


0.75 


100 


50 


0.80613 


0.1907 


0.72854 


8.6001E-1 


9.4777E-1 


700 


50 


0.2179 


0.23669 


0.021002 


1.3836E-1 


3 . 2233E-1 


1300 


50 


0.051488 


0.15441 


0 


9 . 2321E-13 


7 . 8509E-3 


1900 


50 


0.0094727 


0.065348 


0 


0.0000E0 


6 . 1062E-160 



POOLED STANDARD DEVIATION ESTIMATE: 0.17357 



b. Contact Fuzed Warhead 



TABLE 3 

CONTACT FUZED WARHEAD 



PARAMETER VALUES 
POPULATIONS 
WEIGHTS 
X LABEL 
Y LABEL 



(100) ; (700) ;(1300) ;(1900) 
CR14 ;CR13 ; CR12 ; CR11 
0 

RANGE IN METERS 
PROBABILITY OF KILL 



POPULATION 

NUMBER 


NO. OF 
POINTS 


Y MEAN 


Y STD DEV 


0.25 


-PERCENTILES- 

0.5 


0.75 


100 


60 


0.017203 


0.0056864 


0.013458 


0.016796 


1.9509E-2 


700 


60 


0.0092551 


0.0053786 


0.0055805 


0.010485 


1 . 2478E-2 


1300 


60 


0.003111 


0.0055073 


0 


0.000038096 


3 . 3590E-3 


1900 


60 


0.00048311 


0.002269 


0 


0 


1 . 0408E-15 



POOLED STANDARD DEVIATION ESTIMATE: 0.0049179 
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BALLISTIC DISPERSION 



TABLE 4 

BALLISTIC DISPERSION 



PARAMETER VALUES: (0) ; (0 .0005) ; (0 .001) ; (0 . 0015) ; (0 .002) 
POPULATIONS : API ;AP2 ;AP3 ; AP4 ; AP5 

WEIGHTS : 0 

X LABEL : ANGULAR DISPERSION IN RADIANS 

Y LAEEL : PROBABILITY OF KILL 

POPULATION NO. OF PERCENTILES 



NUMBER 


POINTS 


Y MEAN 


Y STD DEV 


0.25 


0.5 




0.75 


0 


100 


0.18605 


0.31466 


0 


2.7032E 


'4 


0.23817 


0.0005 


100 


0.17568 


0.28455 


0 


5.7949E 


"6 


0.27246 


0.001 


100 


0.18497 


0.31487 


0 


6.9311E 


"5 


0.231 


0.0015 


100 


0.1703 


0.2776 


0 


1.9134E 


"5 


0.23033 


0.002 


100 


0.16013 


0.28406 


0 


3.8859E 


"6 


0.13654 



POOLED STANDARD DEVIATION ESTIMATE: 0.29559 



3. RANGE EVALUATION 



TABLE 5 

RANGE EVALUATION 



PARAMETER VALUES 
POPULATIONS 
WEIGHTS 
X LABEL 
Y LABEL 



( 100) ; (700) ; (1300) ; (1900) 
RG14 ;RG13 ;RG12 ;RG11 
0 

RANGE IN METERS 
PROBABILITY OF KILL 



POPULATION NO. OF PERCENTILES 



NUMBER 


POINTS 


Y MEAN 


Y STD DEV 


0.25 


0.5 


0.75 


100 


50 


0.80613 


0.1907 


0.72854 


8.6001E-1 


9 . 4777E-1 


700 


50 


0.2179 


0.23669 


0.021002 


1 . 3836E- 1 


3 . 2233E-1 


1300 


50 


0.046666 


0.15413 


0 


2 . 0463E-9 


2 . 1222E-3 


1900 


50 


0.0094727 


0.065348 


0 


O.OOOOEO 


6 . 1062E-16 



POOLED STANDARD DEVIATION ESTIMATE: 0.16975 
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4. EFFECTS OF TARGET MANEUVERS 



a. Maneuvering Target 



PARAMETER VALUES 
POPULATIONS 
WEIGHTS 
X LABEL 
Y LABEL 

POPULATION NO. OF 



TABLE 6 

MANEUVERING TARGET 



(100) y (700) ; (1300) ; (1900) 
RG14 ;RG13 ;RG12 ;RG11 
0 

RANGE IN METERS 
PROBABILITY OF KILL 



-PERCENTILES- 



UMBER 


POINTS 


Y MEAN 


Y STD DEV 


0.25 


0.5 


0.75 


100 


50 


0.80613 


0.1907 


0.72854 


8 . 6001E-1 


9.4777E-1 


700 


50 


0.2179 


0.23669 


0.021002 


1 .3836E-1 


3 . 2233E-1 


1300 


50 


0.051488 


0.15441 


0 


9 . 2321E-13 


7 . 8509E-3 


1900 


50 


0.0094727 


0.065348 


0 


0.0000E0 


6 . 1062E-16 



POOLED STANDARD DEVIATION ESTIMATE: 0.17357 



b. Non-maneuvering Target 



TABLE 7 

NON-MANEUVERING TARGET 



PARAMETER VALUES 
POPULATIONS 
WEIGHTS 
X LABEL 
Y LABEL 

POPULATION NO. OF 



(100) ; (700) ; (1300) ,- (1900) 
NOM1 ;NOM2 ;NOM3 ;N0M4 
0 

RANGE IN METRES 
PROBABILITY OF KILL 



-PERCENTILES- 



UMBER 


POINTS 


Y MEAN 


Y STD DEV 


0.25 


0.5 


0.75 


100 


50 


0.97438 


0.009994 


9 . 67 51E-1 


0.97609 


0.98113 


700 


50 


0.46539 


0.25853 


2 . 2858E-1 


0.41486 


0.6943 


1300 


50 


0.16126 


0.27431 


6 . 3071E-11 


0.0013258 


0.28531 


1900 


50 


0.012354 


0.075384 


0.0000E0 


0 


0 



POOLED STANDARD DEVIATION ESTIMATE: 0.19231 
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5. EFFECTS OF RADAR ERROR 



TABLE 8 

RADAR ERROR 



PARAMETER VALUES 
POPULATIONS 
WEIGHTS 
X LABEL 
Y LA3EL 



<0) ; (50) ; (100) ; (ISO) 
RD21 ;RD22;RD23;RD24 
0 

PERCENT RADAR ERROR 
PROBABILITY OF KILL 



POPULATION NO. OF PERCENTILES 



NUMBER 


POINTS 


Y MEAN 


Y STD DEV 


0.25 


0.5 


0 


60 


0.3018 


0.41684 


1.4979E-12 


0.014793 


50 


60 


0.2246 


0.33077 


0.0000E0 


0.0000149 


100 


60 


0.13527 


0.25116 


1 . 3878E-17 


0.00014144 


150 


60 


0.083055 


0.1728 


0.0000E0 


0.0017942 



POOLED STANDARD DEVIATION ESTIMATE: 0.29038 



0.75 

0.804 

0.44578 

0.12566 

0.078821 



6. EFFECTS OF MUZZLE VELOCITY AND RATE OF FIRE 
a. Muzzle Velocity 



TABLE 9 

MUZZLE VELOCITY 



PARAMETER VALUES 
POPULATIONS 
WEIGHTS 
X LABEL 
Y LABEL 



(1000) ;( 1500) ; (2000) ;(2500) 
MV1 ;MV2 ;MV3 ;MV4 
0 

MUZZLE VELOCITY(METERS/SEC) 
PROBABILITY OF KILL 



POPULATION NO. OF PERCENTILES 



NUMBER 


POINTS 


Y MEAN 


Y STD DEV 


0.25 


0.5 


0.75 


1000 


60 


0.11223 


0.23417 


0 


0.000029386 


0.055264 


1500 


60 


0.24198 


0.29665 


0.0022309 


0.1129 


0.39962 


2000 


60 


0.35197 


0.27194 


0.09949 


0.31493 


0.57484 


2500 


60 


0.51254 


0.2572 


0.28252 


0.5894 


0.70919 



POOLED STANDARD DEVIATION ESTIMATE: 0.26596 
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b. Rate of Fire 



TABLE 10 

RATE OF FIRE 



PARAMETER VALUES 
POPULATIONS 
WEIGHTS 
X LABEL 
Y LABEL 



(400) ; (600) ; (800) ; (1000) 
FR21 ;FR22 ;FR23 ;FR24 
0 

FIRING RATE (RDS/MIN) 
PROBABILITY OF KILL 



POPULATION NO. OF PERCENTILES 



NUMBER 


POINTS 


Y MEAN 


Y STD DEV 


0.25 


0.5 


0.75 


400 


60 


0.1618 


0.25813 


0.00038014 


0.021753 


0.19313 


600 


60 


0.21132 


0.27047 


0.00079402 


0.058486 


0.39343 


800 


60 


0.30698 


0.30679 


0.0028055 


0.23851 


0.55181 


1000 


60 


0.40055 


0.32489 


0.039921 


0.33058 


0.74919 



POOLED STANDARD DEVIATION ESTIMATE: 0.29132 
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